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Abstract. We study the problem of placing effective upper bounds 
for the number of zeros of solutions of Fuchsian systems on the Rie- 
mann sphere. The principal result is an explicit (non-uniform) upper 
bound, polynomially growing on the frontier of the class of Fuchsian 
systems of a given dimension n having m singular points. As a func- 
tion of n, m, this bound turns out to be double exponential in the 
precise sense explained in the paper. 

As a corollary, we obtain a solution of the so called restricted in- 
finitesimal Hilbert 16th problem, an explicit upper bound for the 
number of isolated zeros of Abelian integrals which is polynomially 
growing as the Hamiltonian tends to the degeneracy locus. This im- 
proves the exponential bounds recently established by A. Glutsyuk 
and Yu. Ilyashcnko. 



1. Zeros of solutions of Fuchsian systems and restricted 
infinitesimal Hilbert 16th problem 

1.1. Fuchsian systems and zeros of their solutions 



Let fl be a meromorphic (rational) n x n-matrix 1-form on the Riemann 
sphere CP 1 with a singular (polar) locus S — {n, . . . , r m } consisting of m 
distinct points. The linear system of PfafHan equations 
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is said to be Fuchsian, if all points Ti E E are first order poles of fi. We 
will denote the collection of all Fuchsian systems of rank n with to singular 
points by & n ,m- Later it will be identified with an m(n 2 + l)-dimensional 
semialgebraic variety. A Fuchsian system can be always expanded as the 
sum of its principal Fuchsian parts, 

(i.2) n = 5>,-, ^- = ^377. '"i^. E^ = ' 

with the residue matrices Ai, . . . , A m (the point r — 00, singular for each 
Ctj , is nonsingular for the sum; if 00 = r rn e Z 1 , then one has to — 1 terms 
corresponding to the finite points n, . . . , t to _i of Z", yet with a nontrivial 
residue A m corresponding to the point at infinity). 

Solutions of a Fuchsian system (1.1) are multivalued holomorphic vector- 
functions on CP 1 \ £, ramified over the singular locus, and growing mod- 
erately (no faster than polynomially in the reciprocal distance to S) when 
approaching points of S along non-spiraling paths. 

We are interested in explicit global upper bounds on the maximal number 
of isolated zeros of linear combinations of the form y(t) — J27=i c i x i(t)i 
t € CP 1 . The initial motivation for this question came from the attempts 
to solve the so called infinitesimal Hilbcrt 16th problem, sec [23, 13] and 
§1.4, but the subject is interesting in itself as well, as linear systems provide 
a vast and essentially unique source of transcendental functions which still 
admit global finite bounds for their zeros [23, §3]. 

Because of the multivaluedness of the solutions, some disambiguation is 
required. Assume that the punctured Riemann sphere CP 1 \ S is covered 
by finitely many simply connected polygonal domains U a (polygonality 
here and below means that each domain U a is bounded by finitely many 
circular arcs and line segments). In each such domain one can consistently 
choose a holomorphic branch of any solution x(t) to (1.1) and form all 
scalar linear combinations y(t.) as above. Clearly, the number 

N a (fl)= sup #{t e U a : c\X\(t) H h c n x n (t) = 0} sC +00, 

x(-), ceC" 

where #{a set} denotes the number of isolated points of the set and the 
suprcmum is taken over all solutions x(-) of the system (1.1), is independent 
of the choice of the branches. 

Let N(Cl, U) be the total number of the zeros counted in this way in all 
domains of a given covering U = {U a }, and the infimum of this total 
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number, taken over all polygonal coverings U, 

N(Sl,U) = ^2N a (n) sc +00, 

(1.3) 

N(f2) = MN(n,U). 

By abuse of language, we say that the number N(f2), if finite, gives a 
global bound for the number of zeros of solutions of the Fuchsian system 
(1.1). Considered as a function 3sf : J^„. m -^NU {+00} on the space of all 
Fuchsian systems (or its subspaces), it will be called the counting function. 

Without additional assumptions the "bound" may well be infinite, 

since even individual linear combinations may have infinitely many isolated 
zeros accumulating to the singular locus S. It turns out that the finiteness 
of the bound is closely related to the spectral properties of the residue 
matrices 

Ai = Rcsfi e Mat„ x „(C), Tj e S, i = l,...,m. 

Ti 

Example 1. — The function y(t) = t 1 + i _I = 2 cos In t is a linear com- 
bination corresponding to the Euler system with U = J) dt in the 
standard affine chart on CP 1 . This combination has infinitely many ze- 
ros accumulating to the two singular points n = and t 2 = 00 with the 
residues A Q = ( _J J) = -A^. 

However, under additional assumptions on the spectra of the residue 
matrices Ai one can guarantee finiteness of the upper bound. The following 
result can easily be obtained by the tools developed in [14, 15]. 

Theorem 1. — Assume that all residue matrices Ai = Res Ti f2 of the 
Fuchsian system (1.1) have only real eigenvalues, 

Spec Ai CM., i = 1, . . . , to. 

Then the corresponding bound for the number of isolated zeros of solutions 
is finite, < 00. 

We will not give the proof of this result since it is purely existential and 
does not allow any explicit estimate for the bound 3sT(0) . Formally it follows 
from Theorem 2 below, whose proof is independent of Theorem 1 . 

1.2. Unboundedness of the counting function 

Denote by y n ,m C ^ n ,m the subspace of the space of Fuchsian systems 
with real spectra of all residues: 

(1-4) J^„. m = {ft e & n ,m - Vr e S, Spec Res ft C E}. 

T 
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By Theorem 1, the counting function N takes finite values on .5f n ,m- Yet 
since ,5f n _ m is non-compact, this finitcncss docs not prevent the counting 
function from being unbounded on S^ n . m . The following examples show 
that this is indeed the case. 

Example 2. — The counting function may grow to infinity together with 
the norms of residue matrices of ft. 

Indeed, for the Euler system, ft = t _1 (^ q) dt with N <G N, the linear 
combination y(t) — t N — 1 has N isolated zeros uniformly distributed on 
the unit circle. In this example iV is the norm of the residues at t = and 
t = oo. 

Somewhat less expected is the fact that the counting function may grow 
to infinity as some of the singular points collide with each other even with- 
out explosion (growth to infinity) of the residues. 

Example 3. - Consider the family of Fuchsian systems il E , e € C, 
£ ^ 0, with the four-point singular locus S = {0, e, 1/e, +00} and four 
residue matrices, each with the real spectrum, 

'0 l\ f \ A l/e = -A, 

0) ' ' W 0) ' A 00 = -A . 



(1.5) 



As the complex parameter e 7^ tends to zero, two pairs of singularities, 
{0,e} and {1/e, 00}, collide. The limit is the Euler system described in 
Example 1 with the residue matrices A +A e and — (^4oo+^4i/e) respectively, 
which has infinitely many isolated zeros accumulating to both singularities. 

This obviously means that limsup| e |^ N(f2 e ) = +00. Indeed, otherwise 
there would exist a finite upper bound N for N(il e ) uniform over all suffi- 
ciently small Consider a triangle K <£ CP \ S which contains at least 
N + 1 isolated zeros of the solution y(t) = t 1 + of the limit system f^o 
strictly inside (cf. with Example 1). By semicontinuity, for all sufficiently 
small values of |e|, the corresponding combination of solutions of the system 
with the matrix fi e will have at least N + l zeros, since the restriction of fi e 
on K depends holomorphically on e. Hence 3Nf(O e ) ^ N+ 1 in contradiction 
with the assumption. 



1.3. Polynomial bounds 

The two examples above show that the counting function ZNT(S~2) grows 
to infinity as approaches the frontier of the subclass y n , m of Fuchsian 
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systems, at least along some parts of this frontier. The main positive re- 
sult of the current work is an explicit upper bound for 3\f(f2) which grows 
polynomially as the Fuchsian system (1.1), identified with its matrix ft, 
approaches the frontier of the class ,5^ n ^ m . 

Here and everywhere in this paper we define the norm of a polynomial 
p G C[xi, . . . , x n ) in several variables as follows, 

(!- 6 ) lbll = ^|ca|, if P{x) = ^c a x Q , a = (a 1 , . . . ,a r ) G Z+. 

The norm \M\ of a matrix M G Mat nX n(C) is taken to be the standard Eu- 
clidean (Hermitian) matrix norm (though any other norm would of course 
be essentially equivalent). 

Let R\, : ,^ n ,m — * be the function defined as follows, 

(1.7) R b (n) = 2+ dist" 1 (t',t")+ |Res T fi|, 

where dist(-, •) is the Fubini-Study distance on CP 1 . The function R\, serves 
to measure the (reciprocal) distance to the frontier of the set & n ,m- Here 
and below R + is a closed subset of nonnegative real numbers. 

Theorem 2. — The counting function N(-) for the number of isolated 
zeros of solutions of Fuchsian systems is explicitly bounded on the subset 
^n,m C & n ,m by a power of the function R\>(-) as follows: 

(1.8) N| V = Vn, m <2° {n ' m) . 

1 ^ n,m 

The constant in 0(n 2 m) is explicit and computable. 

Remark 1. — Note that the Fubini-Study distance on CP 1 is not confor- 
mally invariant, whereas the counting function N obviously is. This may be 
used to improve the bounds: for instance, if m = 3, then any three points 
can be placed at a distance 1 from each other by a suitable conformal 
change of variable t. 

Remark 2. — The term 0(n 2 m) appearing in the bound above is asymp- 
totically equivalent to the dimension of the variety ^n^m- A careful inspec- 
tion of the proof will reveal that this quantity in fact plays the key role in 
determining the magnitude of the upper bound. On suitably defined subsets 
one may apply the same techniques to establish similar bounds depending 
on the dimension of the set, though we don't explore this direction here. 

The proof of Theorem 2 is based on the controlled process of reduction 
of the system (1.1) to a scalar nth order differential equation along the 
lines suggested by A. Grigoriev [8, 9], see also [24]. As a result, we obtain 
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an explicit uniform upper bound for the variation of argument of all linear 
combinations along arcs sufficiently distant from the singular locus of this 
equation. 

By the argument principle, this implies an explicit upper bound for the 
number of isolated zeros away from the singular locus £ of the initial 
system. As for the zeros arbitrarily close to S, we use the method developed 
in [20] for Fuchsian equations. 

1.4. Applications for the infinitesimal Hilbert 16th problem 

The principal motivation for Theorem 2 was the so called restricted infinitesi- 
mal Hilbert sixteenth problem as it was formulated in [6, 11]. 

Let H G WL[x, y] be a real polynomial of degree n + 1 and to = p Ax + q Ay a 
real polynomial 1-form of degree n = max(degp, deg q). The problem is to find 
explicit upper bound for the number Iui,h of real ovals S t C {H — t} on the 
nonsingular affine level curves of H, such that the Abelian integral 

(1.9) Vu,,H(t) = j> u, te (o,6), 

J s t 

vanishes (assuming that such ovals are isolated). The supremum 

(1.10) 3(H) = sup I U , H 

deg uj^n 

depends only on H and is always finite, so that 3 can be considered as a function 
on the linear space J^ n = {deg-ff ^ n + 1} of real polynomials of the specified 
degree. 

In [14, 21] it was shown that the value 3(H) admits a uniform bound, 

(1.11) sup 3(H) = 3„ < +oo, 

yet the proof is non-constructive and the growth of 3 n as a function of n is 
completely uncontrollable. 

All known explicit bounds for 3(H) are non- uniform and finite only for a resid- 
ual (open dense) subset of polynomials H for n > 3. The only case where a 
uniform bound was known to be extendable to a class of degenerate Hamiltoni- 
ans, is that of hyperelliptic integrals, see [18]. Yet the bound in this case is given 
by a tower function (iterated exponent) of unspecified height (iteration depth). 

As a corollary of Theorem 2 we obtain a polynomial global upper bound for 
the number of zeros of Abelian integrals on the same set More precisely, 

consider the expansion of H as the sum of the principal homogeneous parts, 

H = H n+1 +--- + H 1 +H , H k (\x,\y) = \ k H k (x,y). 
Let be a semialgebraic subset in defined by the following conditions, 
(1) the principal homogeneous part H n +i is nonzero: 

H n+1 (x,y)= ^ Cijx'y'^O; 

i+j=n+l 
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(2) the principal part is square-free, i.e., H n+ i factors as a product oi n + 1 
pairwise different linear terms corresponding to different points fi , . . . , 
Cn+i of the projective line CP 1 ; 

(3) the polynomial H itself has exactly n 2 different critical points in C 2 
with pairwise different critical values t\ , . . . , r n 2 £ C (the corresponding 
critical points will automatically be nondegenerate, i.e., the Hessian of 
H is nonvanishing at each of them). 

Obviously, .y? n is a complex affine space and a semialgebraic subset in this 
space. Let : J( n — > R + , H i— » Ri,{H) be the function defined in terms of the 
homogeneous decomposition of H as follows, 

II TTU " + 1 " 2 ™ 2 

(1.12) R^H) = 2 + JLL + J2 dist- 1 ^, 0) + J2 1^1 + E \ T > - ^l" 1 

This is a continuous semialgebraic function on the set jft n with the infinite limit 
on the frontier cL# n . 

Theorem 3. — Tiie number of isolated zeros of Abelian integrals admits a 
computable semialgebraic bound on J% n : 

(1.13) 3(H)\^ n ^Rl v = v n ^2°^\ 

As in Theorem 2, the constant in 0(n 4 ) is explicit and computable. 

Remark 3. — Very recently A. Glutsyuk and Yu. Ilyashenko in a series of 
papers [6, 5, 4, 3] achieved another bound for the number of zeros of Abelian 
integrals, based on completely different ideas. Their bound is finite on the same 
residual subset ^£ n , yet for any fixed n it grows exponentially near some parts 
of the boundary cL# n C ,34f n - 

Theorem 3 is a corollary to Theorem 2 and follows immediately from the fact 
that Abelian integrals satisfy a hypergeometric Picard-Fuchs system of linear 
ordinary differential equations explicitly constructed in [19]. The interested reader 
may easily restore the details of the derivation. 

We do not give a detailed proof of this Theorem here, since in a separate publi- 
cation [2] we establish a uniform explicit bound for the number of zeros of Abelian 
integrals, thus solving the unrestricted infinitesimal Hilbert 16th problem. 
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2. On the general nature of non- uniform bounds 

The form of the bound established in Theorem 2 begs for further analysis, 
as one can question the form in which the numerous different parameters de- 
termining the Fuchsian system (1.1) arc incorporated into a single function 
R\,(£i). In this section we discuss the general form of explicit non-uniform 
bounds on semialgebraic parameter spaces, of which (1.8) is a particular 
case, and study their universality. 

This allows us to conclude that the double exponential dependence of the 
bound (1.8) is independent of the particular form chosen for R\,, as long as 
certain "natural" conditions are satisfied. Furthermore, we show that the 
freedom of choosing the "natural" inverse distance R\, to the frontier will 
affect the powers u n , changing them by a polynomial in n and m. Thus the 
exponential bound for the powers v n ,m (double exponential with respect to 
Rb) cannot be substantially improved in the class of non-uniform bounds 
on the whole class of Fuchsian systems. We explain the precise meaning of 
the "natural bounds" in this section. 

2.1. Semialgebraic sets and their carpeting functions 

Recall that a subset Z C R N of the real affine space is called semialge- 
braic, if it is a finite union of subsets each defined by finitely many polyno- 
mial equalities and inequalities of the form {f a (x) — 0, gp(x) < 0}, where 
fa, gp € R[xi, . . . ,Xn] are real polynomials. A projective subset Z C RP N 
is semialgebraic, if it is semialgebraic in some (hence any) affine chart. A 
subset in is semialgebraic, if it is semialgebraic in the corresponding 
"realification" R 2N ; etc. 

The category of semialgebraic sets is stable by set theoretic operations (fi- 
nite union, intersection, complement) as well as by projections (as asserted 
by the celebrated Tarski-Scidcnberg principle). As a result, we conclude 
that any formula involving quantifiers and polynomial expressions, defines 
a semialgebraic set. For the same reasons the image of a semialgebraic 
set by a polynomial or semialgebraic map is again semialgebraic (a map is 
semialgebraic if its graph is a semialgebraic subset of the Cartesian product 
of the domain and the range) . 

The space J^„. m of Fuchsian systems can be identified in several ways 
with semialgebraic subsets in suitable affine spaces, see §2.2 and Remark 4. 
Following the ideas of A. Khovanskii [15], we introduce the notion of a 
carpeting function for noncompact semialgebraic subsets. Loosely speaking, 
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the carpeting function plays the role of the "reciprocal distance to the 
frontier of the set" . 

Definition 1. A carpeting function of a noncompact semialgebraic 
subset Z is a continuous semialgebraic positive function R: Z — > E + which 
is proper as a map from Z to R + . For technical reasons we will always 
require in addition that i? > 2 on Z. 

In other words, a continuous semialgebraic function R is carpeting for 
a non-compact set Z if R(z) tends to infinity along any sequence {zkj^—i 
without accumulation points in Z ( "converging to a frontier of Z" ) . 

Example 4. — The function R\, is "almost carpeting" on & n , m : it is 
proper, yet not semialgebraic, since the Fubini-Study distance is not alge- 
braic. Yet the latter circumstance is purely technical: one can easily con- 
struct a semialgebraic distance on CP 1 and use it in (1.7) instead of the 
Fubini-Study distance. 

Any two carpeting functions on the same set are related by a two-sided 
Lojasiewicz-type inequality. 

Lemma 1 . — For any two carpeting functions R\ , i?2 on the same semi- 
algebraic set Z, there exist a finite positive constant s such that 

(2.1) VzeZ R 1 2 /s (z)^R 1 (z)^R s 2 (z). 

Proof. — Consider the joint graph {(z, u, v) G R n+2 : z e Z,u = Ri(z), 
v = R 2 (z)}, which is a semialgebraic set by construction. Its projection 
S C M 2 on the (w,u)-plane parallel to the z-direction is semialgebraic by 
the Tarski-Seidenberg principle. 

The function tp(u) = sup{v e R + : (u,v) G S} takes only finite values. 
Indeed, the continuous function R 2 attains its maximum on the compact 
subset {zeZ: = u}. Since the function R\ is carpeting, it assumes 

arbitrarily large values, hence the function ip(u) is defined for all suffi- 
ciently large u € [2, +00). Being semialgebraic by construction, it grows no 
faster than polynomially, <p(u) ^ Cu s , C,s E R+. The constant C can be 
absorbed into the increased power, ip(u) < u s+c , c = log 2 C, for all u > 2. 

The inequality in the other direction is obtained analogously. □ 

In the future any two functions R\ , R2 constrained by the inequalities 
(2.1) on their common domain, will be referred to as polynomially equiva- 
lent. 

Remark 4. — The class of Fuchsian systems (or, more precisely, its 
interior) can be alternatively described as follows. Consider the rational 
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matrix function 
(2.2) 

m— 1 m— 1 

A(t,A)=P(t)/Q(t), P(t)=^M fe t fc , Q{t) = t m +Y,^t\ 

o o 

with the (n x n)-matrix coefficients Mo, . . . , M m _i and scalar coefficients 
Co, ... , c TO _i G C denoted by A e C d , d= (m— l)(n 2 + 1). Assume further 
that P(t) and Q(t) are coprimc. The corresponding linear system will be 
Fuchsian provided that the denominator has no multiple roots, that is, the 
discriminant A e C[A] of the polynomial Q (an explicit polynomial in A) 
does not vanish (the point at infinity is always Fuchsian). 

There is an obvious carpeting function on the space of the parameters, 

(2-3) J?„(A) = 2+^| +||P|| + ||Q||. 

One can easily establish a two-sided inequality between the functions R\, 
and and show that they are polynomially equivalent (one has to express 
the discriminant via differences of the roots of the polynomial Q). Moreover, 
one can easily show that the in upper bound (1.8) the function R\, can be 
replaced by without changing the growth rate of the exponent v, cf. with 
Corollary 1 below. 



2.2. Product structure of the space JP n ,m of Fuchsian systems 

Representation (1.2) together with the special role played by the spectra 
of the residues and geometry of the singular locus, suggests that the space 
■^n,m of Fuchsian systems is naturally (bijectively) parameterized by a 
semialgebraic variety having a rather special product structure: 

<^n,m — ^ m X &n,m — (CP ) X ) , 

(2.4) & m = {(r 1; . . . , T m ) : n G CP 1 , n ± Tj for i ? j} c (CP 1 )" 1 , 
^ m - {{A lt . . .,A m ) : Ai ± 0, Ai G Mat„(C)} C (C" 2 ) m . 

This parametrization is continuous (even biholomorphic) in the sense that 
convergence of the tuples (n, Ai) G ^ m x 3?, n , m in the natural (geometric) 
sense implies the uniform convergence of the corresponding rational matrix 
1-forms on CP 1 on compact sets disjoint with the singular loci of the forms 
(but not necessarily vice versa). 

From now on we will identify the variety ^ n ,m of Fuchsian matrix 1- 
forms on CP 1 having m distinct poles, with points (2m-tuples) from the 
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space (CP 1 )™ x (C™ 2 \ {0})"\ The set y n ,m becomes then a relatively 
closed subset of & n ,m- 

The compactification of the space & n , m preserving the aforementioned 
product structure is the product of projective spaces: 

(2.5) W n , m = (CP 1 )" 1 x (CP" 2 ) m . 

The frontier d-^ n , m = & n ,m \ ^n,m of the variety of Fuchsian systems 
consists of components of three types, 

(2.6) 8^ n 
where: 

(1) the vanishing frontier '1 / n ^ rn consists of the 2m-tuples with one or 
more zero residue matrix A t equal to zero, 

(2) the explosion frontier S n .m is the subset corresponding to one or 
more "infinite" residues A i: and 

(3) c £ m = is the collision frontier corresponding to the union of the 
diagonals \Ji^j{ T i = r j }. 

The parametrization of Fuchsian 1-forms by poles and residues extends 
continuously on the vanishing and the collision frontier, but not on the 
explosion frontier. Note that while the passage to limit as fi tends to S n ^ m is 
impossible, the simultaneous limit on the intersection <^ im nKi,m may well 
make sense in the class of rational systems. One can argue that generically 
such "collision with exploding residues" results in creation of an irregular 
linear system, thus the intersection S n ,m H ^n.m may well be considered as 
an irregularity frontier. 

Examples 2 and 3 suggest that the counting function N: y n ,m — » K+ 
necessarily grows to infinity near the explosion and the collision frontier 
components. On the contrary, Theorem 2 shows that the counting function 
is bounded near the vanishing frontier. 

Together with Theorem 1 this allows to argue that, as long as the problem 
of counting zeros of solutions is considered as an algorithmic mass problem 
formulated for the entire class of general Fuchsian systems, the natural 
domain for the counting function should be the partial closure of the set 
<^n,m, obtained by adjoining the vanishing frontier component y niTn , 

^n, m = {fa, . . . , r m ) : n ^ Tj } x {(A^ A m ) : Spec Ai C K} 

C &*, m = & m x (C" 2 ) m . 

The partial closures m and S / '* m are semialgebraic varieties "parame- 
terizing" the class of Fuchsian systems with to singularities and eventually 



SUBMITTED ARTICLE : POLYFUCHSFINAL.TEX 



12 



GAL BINYAMINI AND SERGEI YAKOVENKO 



zero residue matrices (resp., the class of such systems satisfying the spectral 
condition (1.4)). 

The function R\, defined by the formula (1.7), is clearly a positive function 
on S fi * m which tends to infinity polynomially as ft approaches the fron- 
tier of ^n.m- This function is not semialgebraic because the Fubini-Study 
distance dist: CP 1 x CP 1 — > R + is not semialgebraic, yet this failure is 
purely technical; clearly, there are semialgebraic distance functions on the 
projective line, all of them equivalent to each other; this would transform 
(1.7) into a genuine carpeting function. However, this choice is by no means 
unique. 

Since all carpeting functions on the same semialgebraic set are equivalent, 
Theorem 2 admits the following reformulation. 

Corollary 1. — For any carpeting function R : ,9"*_ m — ► R + the count- 
ing function K(f2) admits an explicit polynomial upper bound of the form 

K(Q)|^. m sC R u (n) for some finite v. 

The constant v depends on n, m and the choice of R. 

2.3. Product spaces and natural carpeting functions on them 

There is no single distinguished (or preferred) carpeting function on a 
semialgebraic variety, and it is even less clear how one may compare carpet- 
ing functions defined on different semialgebraic varieties. Thus it is rather 
difficult to analyze the asymptotic behavior (in n, m) of different bounds 
in invariant terms. 

However, in the particular case of the varieties of Fuchsian systems & n ,m 
one may argue that symmetry considerations distinguish a certain class 
of natural carpeting functions. More precisely, we introduce two special 
classes of product semialgebraic spaces, provisionally called Fermi and Bosc 
types (respectively, with or without the explicit prohibition of coincident 
terms in the product) of several copies of the same semialgebraic variety Z. 
Then we show that on such product spaces one can introduce "anonymous" 
carpeting functions that do not depend explicitly (in the accurate sense 
introduced below) on the number of the copies. These natural carpeting 
functions depend on the product structure types and are defined modulo 
parametric polynomial equivalence (as described in Definition 5). 

Let Z be an arbitrary semialgebraic set. 



ANNALES DE L'INSTITUT FOURIER 



POLYNOMIAL BOUNDS FOR OSCILLATION OF FUCHSIAN SYSTEMS 13 

Definition 2. — A Bose product space Z n is the (Cartesian) power of 
Z, the space of all tuples of points (z\, . . . , z n ), Zi € Z. 

A Fermi product Z* n is the space of all pairwise different tuples of points 
(zi, . .., z n ), Zi e Z: 

Z *n = Z n^\J {z . = Zj} 

Assume now that R : Z — > R + is a carpeting functions on Z and R' 
is a symmetric carpeting function on the Fermi square Z* 2 , i?'(zi,z 2 ) = 
R'(z2,z-\). Starting from these two functions, we can construct carpet- 
ing functions on arbitrary Bose/Fermi products, aggregating these "basic" 
functions by an arbitrary semialgebraic operation. 

Definition 3. — A continuous commutative and associative binary op- 
eration : R + x R + will be called semialgebraic, if its graph is a semialge- 
braic subset in R+, and carpeting, if all sets {(x,y) € R+: x y ^ C} 
are compact and exhaust R + x R + as C — > +oo (i.e., if the function 
f(x, y) = (x y) + 2 is carpeting in the sense of Definition 1 ). 

Because of the commutativity and associativity, the expressions 0™ =1 c, 
make sense for any unordered collection of positive numbers c$. 

Example 5. — The sum x + y, maximum max(x,y), "radius" a/x 2 + y 2 
are semialgebraic carpeting operations. The product is semialgebraic albeit 
not carpeting, yet the "shifted product" x Qy — (2 + x)(2 + y) is. 

More generally, let h: R + — » E + be any semialgebraic homeomorphism 
of M + and any semialgebraic carpeting binary operation (e.g., one from 
the above list). Then the binary operation (x,y) h^ 1 (^h(x) h(y)) is 
again carpeting. 

Definition 4. — A natural family of carpeting functions on the Bosc 
products Z n (resp., Fermi products Z* n ) is a family of carpeting functions 
defined by the formulas 

n 

(2.8) R(z 1 ,...,z n ) = Q)R(z i ), resp., R'(z u ...,z n ) = QR'(zi,Zj), 
j=i i^j 

where R (resp., R') is an arbitrary semialgebraic carpeting function on Z 1 
(resp., on Z* 2 ) and any semialgebraic carpeting commutative associative 
binary operation on the nonncgative ray R + . 

Combining arbitrary natural carpeting functions for the Bose powers 
&n,m — (C n ) m and Fermi powers & m = (CP 1 )*™, by the carpeting 
operation 0, we obtain the class of natural carpeting functions on the 
product spaces ^ n . m — x &n,m as in (2.4). 



SUBMITTED ARTICLE : POLYFUCHSFINAL.TEX 



14 



GAL BINYAMINI AND SERGEI YAKOVENKO 



2.4. Polynomial equivalence of natural carpeting functions 

The degree of freedom used to introduce the natural carpeting func- 
tions on the product spaces is considerable, as one can vary both the basic 
functions R,R' as well as the semialgebraic operation 0. However, these 
variations result in controllable change of the growth rate. 

Definition 5. — Two families of carpeting functions Ri,R 2 : Z n — > R + 
defined on the same family of semialgebraic sets Z n indexed by a natural 
parameter n (in general, the dimensions dim Z n grow to inhnity with n), 
are said to be polynomially equivalent, if there exists a sequence of finite 
positive constants s n G K+, such that, 

(2.9) R\' Sn {z) ^ Ri(z) < R s 2 n (z) VzeZ n ; 

cf. with (2.1), and this sequence grows at most polynomially as n — > oo: 

(2.10) 3c, r < +oo such that s n < cn r , Vn= 1,2,3,... 

An obvious modification allows one to also speak of polynomially equiv- 
alent carpeting functions for families Z a indexed by a multiindex a — 
(ai, . . . , a n ). In this case the respective exponents s a G R+ should grow 
no faster than c \a\ r for some Gnite c, r G K. 

Clearly, there is no reason why two arbitrary families of carpeting func- 
tions be polynomially equivalent: even if this happens, one can involve one 
of the functions into a sufficiently quickly growing sequence of powers and 
destroy the initial equivalence. Yet any two natural parametric families of 
carpeting functions on a family of product spaces turn out to be always 
equivalent in the sense of Definition 5. This can be considered as a para- 
metric version of the Lojasiewicz inequality. 

In complete analogy with Lemma 1, one can show that any two carpeting 
operations are equivalent, in particular, for any such operation 0, 

(2.11) 3 7 < +oo, Va,6 > 2, [ab) lh < a b < (a6) 7 . 

Iterating this inequality, we obtain the following result. 

Lemma 2. — For any carpeting operation as in Dehnition 3, there 
exist two hnite positive constants c, r such that for any n G N and any 
nonnegative numbers xi, . . . , x n > 2 

(2.12) x 1 Q---Qx n < (xi ■ ■ ■ x n ) cnT . 
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Proof. — It is sufficient to prove the inequality for the number of terms 
being a power of two (for the intermediate cases one can replace some of 
the extra terms by units). 

For n = 2™ the proof goes by induction, iterating (2.11). Indeed, if 
7 < 2 r , then the ©-product of 2 m+1 terms is bounded by the product as 
follows, 

(xf mT x? mT y ^ (x 1 x 2 y c2mr ^ (x 1 x 2 ) c2(m+1) \ 

where X\ = X\ ■ ■ ■ x 2 m and X 2 — x 2 m +l ■ ■ ■ x 2 ( m +i) stand for the products 
of the first and the last 2 m terms. The base of the induction can always be 
guaranteed by a sufficiently large c > 0. □ 

Lemma 3. — Any two natural carpeting families on Bose (resp., Fermi) 
powers Z n (resp., Z*n) of the same semialgebraic set Z = Z 1 are polyno- 
mially equivalent to each other. 

Proof. — The proof follows immediately from Lemma 1 and Lemma 2. 

□ 

The function R\, defined by formula (1.7), represents the class of natural 
carpeting functions on the product spaces m and its relatively closed 
subset ^nm- From Lemma 3 we conclude that the assertion of Theorem 2 
does not depend on the specific choice of the expression (1.7); the count- 
ing function N(-) for zeros of solutions of Fuchsian systems admits semi- 
algebraic bound which is double exponential in the complex dimension 
d = dim J?^ m ~ n 2 m of the corresponding variety of Fuchsian systems. 

The following result is the most general form of the polynomial bound 
on zeros of solutions of Fuchsian systems. 

Theorem 4. — For any natural carpeting function R(-) on the class of 
Fuchsian systems considered as a product space m , 

(2.13) N \ y , m < R 2 ° (d \ d = dim c J^ m . 

The constant in 0(d) is explicit, though it depends on the choice of the 
natural carpeting function. 

Proof. — This follows immediately from Theorem 2 since the carpeting 
functions (1.7) form a natural carpeting family, and the polynomially grow- 
ing terms s„ iTO of the polynomial equivalence are obviously absorbed by the 
exponential 2°( m2 "). □ 

Remark 5. — Note that the carpeting function R$ introduced in (2.3), 
is not natural on the space ^ n ,m with respect to the product structure of 
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the latter as it was introduced in Definition 4. Despite that, the function R$ 
is polynomially equivalent to the function R\, in the sense of Definition 5. 

2.5. Near optimality of the double exponential bounds 

The proof of Theorem 4 suggests that the double exponential bound 
asserted in (2.13) is almost optimal. Indeed, a more moderate growth rate 
(say, with the exponent v in (1.8) growing polynomially in n, m, i.e., just 
one notch slower than the exponent) would already be highly dependent on 
the choice of carpeting function. This observation suggests that, at least as 
long as the entire Fuchsian class .y* m is concerned, one cannot expect very 
significant improvements to the (non-uniform) upper bounds of Theorem 4. 

Of course, the situation changes dramatically if instead of the entire 
class ^ n ,m we consider its various subclasses. In [2] we show that for a 
isomonodromic regular quasiunipotent classes of linear systems on CP 1 
one can produce a uniform upper bound for the number of zeros of isolated 
solutions. This bound also turns out to be double exponential, yet because 
of the uniformity one can easily ask about its possible improvement. 

The rest of the paper is devoted to the detailed proof of Theorem 2. 

3. From Fuchsian systems to Fuchsian equations with a 

complexity control 

It is now well understood that there is a substantial difference between 
systems of first order linear ordinary differential equations and (scalar) 
high order linear ODE's with respect to studying zeros of solutions and 
their linear combinations, see [17] and references therein. In this section we 
prove that any linear combination of solutions of a Fuchsian system (1.1) 
from the class J^"* m satisfies a linear equation of order n with explicitly 
bounded coefficients. 

Lemma 4. — An arbitrary linear combination of coordinates of a solu- 
tion of a Fuchsian system with the matrix 1-form S m , in a suitable 
afhne chart t on CP 1 , satisfies a linear ordinary differential equation with 
polynomial coefficients, 

ao(t) y (n) + ■■■ + a n ^(t) y' + a n (t) y = 0, 

a ,...,a n e C[t], gcd(a , • • • ,a„) = 1 

and the following qualifications: 
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(a) deg dj < d — n 2 m, 

(b) \\ aj \\ sC Hooll • R v , with R = R b (fl) and v < 2°(«, 

(c) solutions of the equation (3.1) are holomorphic outside the Unite 
locus £ — {n, . . . , r m } C C such that 

(3.2) |rj| < R, \Ti-Tj\^l, 

in particular, solutions are holomorphic in the disk {\t\ > R}. 

The proof of this Lemma is largely parallel to the exposition from [24, 
§2]. For convenience of the reader we reproduce it with the required modi- 
fications. In what follows we will refer to the ratio 

(3.3) LD = max 



j=l,...,n ||a || 

as the slope of the differential operator I? as in the left hand side (3.1) (and 
the corresponding linear ordinary homogeneous equation). Thus Lemma 4 
asserts an explicit global upper bound on the slope of the equation (3.1) de- 
rived from the Pfaffian system dx = fix, in terms of the carpeting function 
i2b(fi). 



3.1. Preparatory remarks 

We start with choosing a convenient affine chart. From elementary geo- 
metric arguments involving spherical areas, it follows that for any m points 
on the Riemann sphere, there always exists a point at least 0(l/y/m)- 
distant from them in the sense of Fubini-Study distance. We choose the 
affinc chart z G C in which this point corresponds to infinity and normalize 
it so that the Fubini-Study metric is given by the form ■ Such map is 

defined uniquely modulo the Euclidean rotation, and we have an obvious 
inequality dist(zi,z 2 ) < \z\ — z 2 |. 

By construction, in the chosen chart z the affine coordinates n , . . . , r m of 
the singular points of a system (1.1) with R\,(Cl) = R, satisfy the inequality 
\ z j\ ^ 0(^/m) for all j — 1, . . . , m. Since all terms in the expression (1.7) 
are nonnegative and the residues are invariant, we conclude that in this 
chart the inequalities 

(3.4) \ Zi \^0{Vm), j = l,...,m, \ Zi - Zj \^l/R, i^j 

hold simultaneously By a suitable affine rescaling t = cz we can normalize 
the singular locus of the system to satisfy the inequalities (3.2). In the 
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affine chart thus constructed, the Pfafhan linear system (1.1) is equivalent 
to a system of linear ordinary differential equations of the form 



Finally, performing a linear transformation x^Ux with a unitary (con- 
stant) matrix U, we may without loss of generality assume that the linear 
combination in question is just the first component, y(t) — xi(t). This 
transformation does not affect the norms of the residues, hence the in- 
equalities (3.4) are preserved. 

3.2. Objects defined over Q and their complexity 

If we consider the entries of the residue matrices Ai and the affine co- 
ordinates Ti G C of the singular points as coordinates on m , denoting 
a general point of this space by A, then the coefficients of the linear sys- 
tem (3.5) become integer numbers, more precisely, only and ±1. In other 
words, the entries of the matrix 1-form Q\ are defined over the field Q on 
the space with coordinates (t, A). In sharp contrast with the fields of real 
or complex numbers, for objects defined over the field Q one can define 
their complexity (bitsize, or, more precisely, the exponent of the bitsize). 
This definition may be introduced by many equivalent ways. For instance, 
as follows: 

(1) Complexity of a polynomial p € Z[A, t] is its norm ||p|| e Z+, the 
sum of all absolute values of its integer coefficients; 

(2) Complexity of a rational function R <G Q(A, t) is the minimal value 
of sum ||p|| + over all representations R = p/q with p,q £ Z[A, t]; 
in particular, complexity of an irreducible rational number r = p/q, 
p,qeZ, gcd(p,q) = 1, is \p\ + \q\; 

(3) Complexity of any object constructible over Q(\,t) is the sum of 
complexities of all rational functions explicitly occurring in the con- 
struction. For instance, complexity of a rational matrix 1-form £1 is 
the sum of complexities of all its entries; complexity of a differential 
operator with rational coefficients is the sum of complexities of its 
coefficients etc. 

This definition (as most constructions from the complexity theory) in fact 
applies not to the object in question, but rather to a specific representation 



(3.5) 



dx 

~d7 




\Aj\^R, \tj\^R, 
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of this object. It will be important in this paper that the objects under 
discussion are described by their explicit representation. 

Whenever two representations are related by a simple transformation 
(e.g., expansion or reduction to common denominator), their complexity is 
usually comparable. More generally, any explicit construction over Q ad- 
mits explicit control over the complexity growth of all intermediate objects. 
This observation is the primary reason for introducing the objects over Q 
in the context where the initial field is C or R. We will apply this technique 
to the straightforward process of derivation of a scalar differential equation 
of high order from a system of first order equations (consecutive derivations 
and elimination). 



3.3. Derivation of the scalar Fuchsian equation over Q(X,t) 

By direct inspection of the formula (3.5) one can see that the matrix 
product P = n™(^ — T i)A(t) is a polynomial matrix with entries in Z[A, t] 
of degree ^ m and norm not exceeding m ■ 2 m ~ 1 ^ 2 m (here and below 
we use the fact that the norm is multiplicative, \\pq\\ ^ • ||<?|| for any 
two polynomials p, q). For brevity, we call elements of the ring Z[X,t] the 
lattice polynomials. 

Lemma 5. — The first component of any solution of the system (1.1) 
considered as a linear system with coefficients from the field Q(A, t), satisfies 
a linear nth order differential equation with lattice polynomial coefficients 
of the form 

oo(A, t) + oi (A, t) t/"- 1 ) + • • • + a n (X, t) y = 0, 

(3.6) 

a , . . . , a n e Z[A, t], a ^ 0. 

The coefficients aj are lattice polynomials in the variables A, t having 
degrees ^ n 2 m and explicitly bounded norms, \\aj\\ s; 2°(™ m \ 

Note that the norms in the formulation of Lemma 5 arc computed with 
respect to the ring Z[A, t], hence they are nonnegative integer numbers. 

Proof of the Lemma. — By virtue of the system (3.5), the higher deriva- 
tives of the function y(t) can be expressed inductively as linear combina- 
tions of the dependent variables x — (x\, . . . , x n ), 

(3.7) V w (t) = &-x, &e(Q(A,i)) m , 

with the n-vectors over the field Q(A, t) determined by the recursive rule 

Co = (1,0,..., 0), 

(3.8) ^eMat„(Q(A,t)), d t = d/dt 
t,k+l = Ct^fc + Sfc • A, 
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(we treat a: as a column vector and £ k as row vectors from the dual space). 

The structure and complexity of the n-covectors can easily be estab- 
lished by the inspection of the formulas (3.8). Denote 

m 

Q(\,t) = H(t-n)eZ{\,t}, ||Q|| = 2"\ dcgQ-m, 

i=l 

the discriminant, so that A — P/Q in (3.8). 

The n- vectors £o> ■ • • > £n defined by (3.8), satisfy the following properties: 

(1) Their multiples a k = Q k (, k arc lattice n-vector polynomials for 
all k = 0,...,n, of degree no greater than km, and the norms 
IKH sC (3nm)"2" m . 

(2) The wedge product of any n vectors from this collection, the poly- 
nomial n-form 

Q n(n+1)/2 £oA---A&_iA& + iA---A£ n) 

has as its only coefficient a lattice polynomial of degree at most 
n 2 m and norm bounded by n!(3nm)" 2 2n m . 

(3) The wedge product £o A • • • A £„_i does not vanish identically on 

,9* 

Indeed, the polynomial 1-forms a k — Q k S,k satisfy the recurrence rela- 
tions 

(3.9) a k+1 = Q ■ d t a k - (kd t Q - P) • a k , 

which follow directly from (3.8). The bound for the degrees dega k is ob- 
vious since degQ,degP < m. The bound for the norms ||afe|| can also be 
derived from the recursive equations, 

K +1 || = \\Q-d t a k -(kd t Q--p)a k \\ 

sC IIQH (dega k ) \\a k \\ +fc(dcgQ) ||Q|| ||a fe || + n||P|| ||a fe || 

s$ (2 m km + 2 m km + 2 m nm) \\a k \\ 

sC 3 • 2 m ■ nm\\a k \\ , 

so that an easy induction gives \\a k \\ ^ (3nm) n 2 nm . To establish the bound 
for the wedge products, we expand the nxn determinants using the Laplace 
expansion involving n! summands, each having the degree ^ \n(n — l)m 
and norm s$ (3nm)" 2 2" 2m 2°^ 2m \ 

Finally, to establish (3), we note that for specific Fuchsian systems (1.1) 
and for specific solutions y(t) we can apriori guarantee that the functions 
y(°\t) , . . . , y^^ 1 ^ (t) are linearly independent. For instance, the function 
y(t) = l+t+- • ■+t n ~ 1 is obviously linearly independent with its derivatives 
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of all orders < n — 1. This function is a linear combination of solutions of 
an Euler system with two singularities at the origin and at infinity, which 
is a particular case of a Fuchsian system. For the corresponding value of 
the parameters A, the Wronskian £ A • • • A £„_i is not identically zero as a 
function of t. 

To derive the equation (3.6), it remains to note that any n + 1 vectors in 
the n-space (Q(A, t))" are linearly dependent over this field. In our partic- 
ular case the vectors £ > ■ ■ ■ , £n-i are linearly independent, so the vector 
expands as a linear combination, = co£o + ' ■ ■ + c ra _i^ n _i, of the previous 
vectors with coefficients Cj <E Q(A, t). These coefficients can be found using 
the Cramer rule, as ratios of the determinants, 

(o -, ^ M-A^i-i A&+i A---A£ n . n 1 

(3.10) c 4 = — — , i = 0,...,n-l. 

4o A • • • A 4„_i 

The numerator and the denominator of any such fraction are rational func- 
tions defined over Q with the denominators being known powers of the 
discriminant Q. Multiplying the linear identity by the common denomi- 
nator, we obtain a linear identity between the n- vectors £o>---i£n with 
lattice polynomial coefficients, which by (3.7) produces a linear nth order 
differential equation for the function y. 

The leading coefficient ao of this equation is a nonzero lattice polynomial, 
therefore its norm is a positive integer number which is necessarily ^ 1 . □ 

3.4. Explicit bounds on coefficients of the scalar equation over 

C(t) 

Specializing the equation (3.6) for each particular value of the parame- 
ters A would produce an equation of the type (3.1) with coefficients from 
C[t]. The bound for the degrees of these coefficients will obviously be the 
same as for (3.6). However, the condition ||&o|lz[At] = 1 does not prevent 
the specialization ao(A, •) of the leading coefficient from having an arbi- 
trarily small norm ||a (A, -)llc[t] or simply vanishing. In this subsection we 
prove, following [9, 24], that this vanishing may happen if and only if all 
other (nonlcading) coefficients also vanish. Moreover, we show that the ra- 
tios ||aj(A, -)llc[t] / || a o(A, -)llc[t] are explicitly bounded as functions of A. 
(Lemma 6 below). 

Consider a linear parametric differential equation (3.6) with lattice poly- 
nomial coefficients ao,...,a n G Z[t, A]. Denote by A C C the singular 
locus corresponding to the parameters A for which the leading coefficient 
ao (A, •) vanishes identically as a univariate polynomial. 
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Solutions of the equation (3.6) may be considered as holomorphic func- 
tions of t, A by fixing an initial condition at some point t = to and treating 
A as parameters of the equation varying near some value Ao, provided that 
the bidisk {\t — t | < a , I A — A | < b} is free from singular points of the 
equation. Yet in general these solutions will exhibit singularities on the lo- 
cus {(X,t): a (\,t) = 0} C C x C which contains some exceptional lines 
{A*} x C for A* e A, while intersecting the other (generic) lines {A} x C, 
A ^ A, only by isolated points. 

Definition 6. — The family of equations (3.6) has an apparent singu- 
larity at a point A* e A C C n , if there exist: 

(1) a small neighborhood U C of A*, 

(2) an open disk D C C, 

(3) a fundamental system of solutions y\ (A, t) , . . . , y n (A, t) of the family, 

such that each function yi(X, •) analytically depends on A in the disk D as 
long as A £ U \ A, and remains bounded in D as A — > A*. 

By the removable singularity theorem, this means that a fundamental 
system of solutions admits uniform analytic limit in D as A — > A*, though 
the linear independence of the functions j/i(A*, •) may fail at the limit. 

Denote by D\ the linear differential operator corresponding to the linear 
homogeneous equation (3.6) after specializing the value A ^ A. Denote by 

(3.11) <p(\) = ZD x = max jj^ "j 11 A e C x A, 

j=i,...,n ||a (A, -)|| 

the slope of this operator. Without loss of generality we may assume that 
the polynomials oo(A, ■),-■■, a„(A, •) arc mutually prime in the ring C[t}: 
nontrivial common divisor may occur only on a proper algebraic subset of 
C which can be combined with A. 

Lemma 6. — Assume that the family of equations (3.6) with lattice 
polynomial coefficients of known degree ^ d and bounded norms \\cij\\ ^ 
M, depending on the complex ^-dimensional parameter A e C, has only 
apparent singularities at all points of the locus A which is a proper algebraic 
subset ofC. Then 

(3.12) VA such that |A| < M, p(A) < M d ° W 

with an explicit constant in O(-). 

Sketch of the proof. — The proof of the Lemma almost literally repro- 
duces that of Lemma 5 from [24] with minimal modifications. 
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More specifically, we start with the observation that the slope ip is locally 
bounded everywhere on <C V \ A. This follows from the assumption that the 
family (3.6) has only apparent singularities, in the same way as in the proof 
of Lemma 4 from [24] . The only required change is to replace the unit disk 
by the disk D from the Definition 6. 

Next, we consider the positive subgraph X of the function ip(-) over the 
ball of radius M. This subgraph is a semialgebraic subset in X C C +1 , 
defined by lattice polynomial equalities and inequalities. Since the function 
ip is locally bounded near each point of the compact ball {|A| ^ M}, the 
semialgebraic set X is globally bounded. 

The lattice polynomials defining X are of degrees not exceeding 0(d) 
in Re A and Im A. In the same way the norms of these polynomials are no 
greater than O(M). Finally, the dimension of the space in which X resides, 
is 0(v). By the quantitative version of the Tarski-Seidenbcrg quantifier 
elimination theorem [7, 10] and, most recently, [1, Theorem 7.2], the diam- 
eter of X does not exceed M d ° { " ] ', as asserted in (3.12). □ 

Proof of Lemma 4. — The proof directly follows from Lemmas 5 and 6: 
we first derive the scalar equation (3.6) over Q(A, t) and transform it to the 
form with lattice polynomial coefficients by eliminating the denominators. 

The resulting family of scalar equations (3.6) has only apparent singular- 
ities, since for any parameter A solutions of the initial system (3.5) depend 
analytically on A in any disk free from singular points S = E\ C C of 
the latter. By Lemma 6, we conclude with the inequality asserted by the 
Lemma, for all values A outside from the vanishing locus A = {X: a (X, ■) = 
0}. This locus is a proper algebraic hypersurface of the complex affine space, 
hence by continuity the inequality ||oj|| ^ ||a || • R v holds everywhere on 
the ball {|A| < R}. □ 



4. Counting zeros of solutions of Fuchsian equations 

In this section we complete the proof of the main Theorem 2. It is based 
on a complex counterpart of the de la Vallee Poussin theorem [22] , which 
asserts that the variation of argument of any (complex) solution of a monic 
linear differential equation (with the leading coefficient identically equal to 
1) along any rectilinear segment can be explicitly bounded from above in 
terms of the magnitude of the non-principal coefficients of this equation. 
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4.1. Lower bounds of polynomials away from their zeros 

We start with an observation which relates the norm (1.6) on univariate 
polynomials with the supremum-norm on disks of finite radius centered at 
the origin. Clearly, for any finite d the two norms on the finite-dimensional 
space of polynomials of degree ^ d are equivalent. However, we shall also 
be interested in the explicit constants involved in this equivalence. 

The first side of this equivalence is easy to describe explicitly. Indeed, 
\\p\\ = 1 implies an easy upper bound \p(t)\ < R d for \t\ < R. We will prove 
the inequality in the opposite direction, establishing a lower bound for \p(t)\ 
for all t € C distant from the null locus of p, provided that ||p||c[t] = 1- 

Lemma 7. — Let p = ^2i PjV be a polynomial of degree d and unit 
norm ||f>||qt] = L an d £={ieC: p(t) = 0} the zero locus of p. 

Then for an arbitrary point t £ S in the disk of radius R > 2, we have 
an inequality 

(4.1) \p(t)\ ^ 2-°^(r/R) d , r = dist(t, S), \t\ sC R. 

Proof. — For a polynomial p = X)j=oPi^ °f um ^ norm an d degree d, at 
least one of its d+ 1 coefficients must be at least l/(d+ 1) in the absolute 
value. Writing the Cauchy formula for the corresponding derivative at the 
origin, we conclude that 

PJ=^tJ r>(C)dC, hence i- / |p(C)| ds^. 
2 ™J\t\=i 2 ™ J\t\=i "+ 1 

By the mean value theorem there exists a point t* on the unit circle, such 
that \p{U)\ > 

Let K C C be the disk of radius 2R, centered at the origin. The polyno- 
mial p can be written as the product, 

(4.2) p(t)=a J] (t-W II (!-£)> aeC > 

\tj\>2R \tj\^2R 

where tj are roots of the polynomial, counted with their multiplicity. Denote 
by d± the number of roots of S inside (resp., outside) the disk {\t\ < 2R}, 
counted with multiplicities so that d_ + d + = d. 

A lower bound for \a\ follows from the inequality |p(t*)| = l/(d+ 1). 
For t = t* the first product in (4.2) is majorized by (2R + l) d - , while the 
second product does not exceed (4/3) d+ , since |t|/|ij| ^ l/(2R— 1) 1/3 
for all terms in it. Altogether we conclude that 

\q\>J-. 

1 1 " d+1 (2R+l) d -~ 
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This implies the lower bound for \p(t)\ with \t\ < R: 
\p(t)\>\a\r d -(l/2) d +. 
Combining these two estimates, we conclude that 

for all 2 s$ R < +00 and R > r > 0. □ 

Remark 6. — One can give a different bound using the Cartan lemma 
[16]. This lemma establishes an explicit lower bound for a monic polynomial 
outside a union of disks of given total diameter. For our purposes it is less 
convenient, since the polynomial in question is not monic (and can actually 
be very far from it). 

4.2. Variation of argument of solutions along arcs 

For the next step, we give an upper bound for the variation of argument 
of complex-valued solutions of a Fuchsian equation along a circular arc (or 
line segment) in terms of the slope of the equation and the normalized 
length of the arc. 

Lemma 8. — Let D be a differential operator of order k with coefficients 
of degree < d and slope S = ZD, and 7 a closed circular arc or line segment 
disjoint with the singular locus S, which belongs to the disk of radius R 
centered at the origin. 

Then the variation of argument of any nonzero solution of the homoge- 
neous equation Dy = along the arc 7 is explicitly bounded, 

(4.3) VarArgy(t)| 7 < kSL(R/r) 0( - d \ 

where L is the length of the arc 7 and r = dist(7, S). 

Proof. — Without loss of generality we may assume that the polynomial 
coefficients of the operator D are normalized as follows, 

||ao||=l, IKKS, j = l,...,k, 

where S = ZD. 

On the disk of radius R the absolute values of the non-principal coef- 
ficients are bounded, |aj(t)| ^ SR d . Restricting the equation Dy = on 
the arc 7 parameterized by the arc-length, we obtain a homogeneous linear 
ordinary differential equation. Dividing by a | 7 and applying Lemma 7 to 
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this polynomial, this equation can be brought into the monic form with 
explicitly bounded complex valued coefficients, 

y W +A 1 (s)y^ + --- + A k (s)=0, se[0,L], 
|A,(s)| s? 2°^ • SR 2d /r d < S{R/r)°^ d \ j = 1, . . . , k. 

Applying the complex generalization of the de la Vallee Poussin nonoscil- 
lation theorem [22, Corollary 2.7] to (4.4), we conclude that the variation 
of argument admits the bound 

(4.5) Var Argy(-)L < Oil) ■ kLA, A= max max|A,-(s)|. 

j=l,...,k s£*y 

Substituting the bounds (4.4) into this inequality, we conclude with the 
bound (4.3). □ 

The assertion of Lemma 8 allows one to place an explicit upper bound 
for the number of isolated zeros of any solution to a linear equation with 
rational coefficients in a simply connected polygonal domain away from 
the singular locus. Indeed, one can estimate the variation of argument of 
any solution along each arc bounding the polygonal domain, and apply the 
argument principle. 

In addition to the explicit bound on the variation of argument along 
arcs distant from the singular points, one can slightly modify the argument 
above to give a bound on the variation of argument along a small circle 
centered at a Fuchsian singular point as well. 

Lemma 9. — Let to G £ be a Fuchsian singular point of a linear oper- 
ator with polynomial coefficients of order k, degree d and slope S, which 
belongs to the disk of radius R > 2 and is at least r-distant from all other 
points t\, . . . , td € E. 

Then the variation of argument of any solution y(t) along the circle 
1p = {I* ~ I = p} is bounded, 

(4.6) VarAr g2/ (-)| 7p ^ kS(R/r)°^ 

uniformly over all sufficiently small p > 0. 

Proof. — Note that the parallel translation p(-) (Tr#p)(-) = p( ■ + R) 
is a linear operator of explicitly bounded norm on the space of polynomials 
of degree < d. Indeed, by the Taylor formula, 

II Tr^|| < £ < |b|| • £ {R f < 2°^R d • |b||. 

i=0 ' 
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Thus without loss of generality we may assume that the Fuchsian singular- 
ity to = is at the origin, and the slope of the new differential equation is 
at most 2°^SR d , where S is the original slope. 

Since the point is Fuchsian, the equation can be written in the form 

(4.7) t k q Q (t) yW(t) + i fc -V(*) V (k - 1] {t) + ■■■ + q k (t)y(t) = 0, 

where qo , qi , . . . , q k are polynomials of degrees ^ d — k,d — k+1, ...,d 
respectively, see [12, Proposition 19.19]. Since the norm is multiplicative, 
and \\t j \\ = 1 we conclude that ||<7j||/|M| < 2°^SR d (see above). By 
Lemma 7, \qo(t)\ ^ 2~°( d V d and, obviously, \qj(t)\ < ||gj|| for all sufficiently 
small \t\ = p. 

Letting t = ps and noting that = P~ k ~^k, we rewrite the equation 
above in the form 

(4.8) s k q (ps) yW(s) + s^q^ps) y^is) + ■■■ + q k (ps)y(s) = 0, 

with \s\ = 1. Notice that the rescaling of precisely cancels out with the 
rescaling of the vanishing terms t k in the coefficients. This circumstance is 
unique to Fuchsian differential operators. 

After dividing the equation above by the leading term (nonzero and 
bounded from below) , we obtain a monic linear ordinary differential equa- 
tion on the unit circle \s\ = 1 with coefficients not exceeding an explicit 
constant A = 2°^ S(R/r) d . The bound for the variation of argument along 
the circle again follows from (4.5). □ 



4.3. Counting zeros in topological annuli 

In this section we recall an "argument principle" for multivalued func- 
tions. This principle asserts that the number of zeros of a solution of a 
Fuchsian differential operator in an annulus is explicitly bounded, provided 
that the monodromy of the operator has all eigenvalues on the unit circle 
and the annulus is at a positive distance from the singular locus. The proof 
is based on an idea due to Petrov (the so called "Petrov's trick"). The 
result we need appears (in slightly less general form) in [20] . 

Let K be an annulus bounded by two disjoint circles C\_i with centers on 
the real axis, and D a monic linear operator with coefficients holomorphic 
in the closure K and real on the real axis K n R (we will refer to such 
operators as real). In this case Lemma 8 asserts that there exists an upper 
bound B = B(D, K) for the variation of argument of any solution of the 
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homogeneous equation Dy = in K along each of the boundary arcs of 
dK, 

(4.9) VarArgyO)!^,, ^B = B(D,K). 

Lemma 10. — If the monodromy of a real differential operator D along 
the equator of an annulus K symmetric around the real axis has all eigen- 
values on the unit circle, then the number of zeros of any solution in K is 
explicitly bounded, 

(4.10) X(D,K) < (2fc + l)(2B + l), 

where k is the order of the operator and B = B(D, K) the bound from 
(4.9). 

Proof. — This Lemma almost literally coincides with Theorem 2 from 
[20]. To simplify the exposition, we first make a real conformal automor- 
phism of the complex line CP 1 , which transforms the given annulus K into 
an annulus K' bounded by two concentric circles centered at the origin, 
and then pass to the logarithmic chart £ on the universal covering over K' . 
As a result, we arrive at the following problem. 

Let II be the "vertical" strip c_ ^ Re ( ^ c + and D a linear ordinary 
operator with coefficients holomorphic in II, 27ri-periodic in this strip and 
real on II n M. Let B be the upper bound for the variation of argument 
of any solution of the equation Dy = along the segment of length 2ir on 
any of the lines Re £ = c± . Because of the periodicity, the shift operator 
M : y(-) i— > y( ■ + 2m) is an automorphism of the space of solutions of the 
equation, corresponding to the monodromy of the original equation in the 
annulus. 

We need to show that if all eigenvalues of M are of unit modulus, = 1, 
then the number of zeros of any solution in the rectangle II — II n {| Im £| ^ 
2n} is at most (k + l)(2B + 1); this is equivalent to the assertion of the 
Lemma. 

By the argument principle, it is sufficient to place an upper bound on 
the variation of argument of any solution / of the equation Df = along 
the perimeter of the rectangle, which consists of two vertical segments on 
the lines Re ( = c± and two horizontal segments 

S± = {Im C = ±2tt, c_ s$ Re C sC c+}. 

The contribution of the vertical segments is explicitly given by the as- 
sumptions of the lemma; it does not exceed 2B. It remains to estimate the 
variation of arguments along the horizontal segments. Each increment of 
Arg/(C) by an angle it or more on a connected interval of ( implies, by 
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an intermediate value theorem, that at some point ( on this interval the 
imaginary part Im(/x/(£o)) vanishes, for an arbitrary choice of the complex 
number ji ^ 0. This observation implies the inequality 

VarArg/| s± < tt(i/± + 1), v ± = {C G S± : Im( M *7(C)) = 0} 

valid for any complex number /x =/= 0. In general, one has no explicit control 
over the imaginary part of a solution along the segments S±. However, if 
/ is real on the real segment Sq = II n R, then by the Schwarz symmetry 
principle, f(( + 2iri) = /(C — 27ri) for ( e S . Therefore for any fi with 

ImI = i, 

Im(^V)|s ± =Im( M ^ 1 M ±1 /)| So =P M /| So , 
where is the Petrov operator 

(4.11) P M = I^M-mM" 1 ). 

This yields an inequality relating the number N(f,H) of complex zeros 
of a real solution / in the rectangle II to the number of real zeros of the 
function /' = P M / on So C II: for any fi with \fi\ = 1, 

(4.12) N(f, n) < i- [2B + 27r(7V(/', So)) + l) < (25 + 1) + iV(.f ,n). 

This inequality is a complex analog for the Rollc inequality relating the 
number of zeros of a real smooth function and its derivative. The inequality 
remains valid also in the case when /' = P M / = 0, if we set -/V(0, II) = 0. 
This follows from the symmetry of the function / (the contributions of 
VarArg/ along the segments S± cancel each other). 

To complete the proof of the Lemma, we show that the product (com- 
position) of operators = P™. ; (where /ij range over the eigenvalues 
of the monodromy operator M and m, over the corresponding multiplici- 
ties) vanishes on the linear space of solutions of the operator Dy = in II. 
This assertion is an analog of the Hamilton-Cayley theorem for the Petrov 
operators. 

Indeed, suppose M is a Jordan block of size m with an eigenvalue \i and 
let C V\ C • • • C V iii C denote the corresponding flag. Then M acts on 
Vj+i/Vj as multiplication by /x, and P^ acts by zero. Thus P^Vj+i C Vj 
and it follows that P™ = 0. It remains to note that Petrov operators with 
different numbers /i commute with each other. 

Now the proof of the Lemma for functions real on K is straightforward: 
if Df — 0, then P_d/ = 0, and by (4.12), the number of isolated zeros of / 
does not exceed k(2B + 1), where k is the order of D, equal to the number 
of the Petrov operators forming the product P. 
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Finally consider the case where / is a general complex valued solution. By 
the same argument as above, the total number of zeros of / in Lin { | Im (\ ^ 
2ir} does not exceed 

(4.13) B+ ± (VarArg/| s+ +VarArg/| s _) . 

Assume that Im/|s ± ^ (if this assumption fails, consider Re/|s ± in- 
stead). Let f± denote the function which agrees with lm('M ± f) on So and 
is extended analytically onto II. The values of on So correspond to the 
imaginary part of / on S* respectively. 

Since D is a real operator, / ± is again a solution of the equation Dy = 0. 
These functions are real on So, the previous considerations apply, and we 
conclude that 

VarArg/| s± ^2nN(f ± ,S ) <fc(2B + l). 
Combining this with (4.13) we finally obtain the bound (4.10). □ 

4.4. Symmetrization 

In order to apply Lemma 10, to a given annulus, one needs the as- 
sumption that the coefficients of the operator D are real on the real axis, 
eventually after the rotation of the axis. In this section we show that this 
additional assumption can in fact always be achieved, if the operator D is 
obtained from a system of linear differential equations (1.1). This construc- 
tion will be referred to as the symmetrization of the operator relative to a 
given axis. 

Let / £ &{U) be a function holomorphic in an open domain U C C. The 
function f\ defined in the domain W by the formula 

(4.14) /*(*) = /(*)> t£U^ = {i:t£U}, 

(the bar stands for the complex conjugacy) will be called a reflection of / 
in the real axis. Obviously, p is holomorphic, p £ &(U^). For polynomials 
the reflection consists replacing all (constant) coefficients by their complex 
conjugates: 

m m 

p(t) = ^c k t k <^=> p f (i) = c k t k , c , ... , c m £ C, 

fc=0 fc=0 

and it extends to the rational functions in the obvious way. 

Let U = W C C be a domain symmetric with respect to the real axis 
RcC. A C-linear subspace L C t?(U) of functions defined in U, is called 
symmetric (with respect to the real line), if together with each function 
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/ £ 0(U) it contains the function p. A symmetric subspace is closed 
under taking the real/imaginary part on LTlR: for any f £ L the functions 
Re f\uns. and Im funs, can be analytically continued from the subset of real 
points of U onto the entire domain U using the formulas Re / = \{f{t) + 
P(t)), Im / = ^(f(t) + P(t)) respectively. A symmetric linear space (in 
case it is finite-dimensional) always admits a basis of functions real on 
U fl K. A monic linear ordinary differential operator D £ &(U)\d\ defines 
a symmetric space of solutions, if and only if its coefficients are real on R, 
i.e., D £ @-r(U fl R)[<9]. Such operators will be referred to as real (on R). 

Any finite-dimensional subspace L C ff(U) of functions can be embed- 
ded into the symmetric completion LP — L + by simply adjoining the 
reflections /' of all elements / £ L. As a consequence, each linear opera- 
tor D £ &(U)[d] can be considered as a "right factor" of a real operator 
D e £ ^r({ ; fl M.)[d] of higher order with real coefficients, which vanishes 
on all functions from L e . However, the slope Z D e of the operator may 
increase uncontrollably (relatively to Z D) after such symmetrization. 

Example 6. — The function f(t) = i + y^ 100 satisfies a first order 
equation 

t" 

y' - i+ wioo y = 

1 ^ 100 1 

which has very small (in the absolute value) coefficient in the disk \t\ < | 
and the slope ZZ? very close to 1. On the other hand, the symmetrization 
L e is generated by the functions 1 and t 100 /100 and is defined by the 
second order equation 

y"-y!/' = 0. 

This equation has uniformly large coefficients in the same disk (and a pole 
at the origin) . Even the slope of the corresponding second order operator 
is relatively large, Z D e = 99. 

However, at the level of linear first order systems the symmetrization 
can be achieved very easily. Let <E ^ n ,m be the rational matrix 1-form 
built from the 1-form £ ^ n ,m on C by the reflection in the real axis 
as above (replacing all coefficients of the entries oJij by their conjugates). 
Then each solution of the reflected system dX^ = • can be obtained 
by reflection of the corresponding solution of the initial system da; = fix 
(1.1) and vice versa. Therefore the symmetric linear space LP spanned by 
all components of a solution of the system 

(4.15) dy = n e y, fi e = , y £ C 2 ™, 
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includes among its solutions the symmetrization of the solution of the ini- 
tial system. Reducing the doubled system (4.15) to a scalar equation, one 
obtains an explicit expression for the operator D°. 

Clearly, one can replace the real axis R C Cby any (real one-dimensional) 
line i in the complex plane C. Any such line can be mapped to the real 
axis by a suitable affine transformation ip: C — > C (in fact, rotation and 
translation suffice). The corresponding reflection operator replacing (4.14), 
can be explicitly written as follows, 

/H/| = (/o V )to^. 

Let now U be an arbitrary concentric annulus free from singular points of 
a Fuchsian system fi <G ^ n ,m- Among all systems obtained by reflection of 
the matrix 1-form fi in different axes I passing through the center, one can 
always find a system (denote it again by fit) whose singular points (except 
for the singular point at the center) will be sufficiently far away from the 
singularities of fi. More precisely, we can always guarantee that the matrix 
form fi e of the symmetrized system (4.15) satisfy the inequality 

(4.16) i?b(fi G ) «S R\,W, v = 0{n + m) 1 

with an explicit constant in O(-). Indeed, the norms of the residues of the 
symmetrized system are majorized by the norms of residues of fi (residues of 
the matrix form fit have the same norm as those of fi) . As for the distances 
between the singular locus £ = {tj} and its mirror image , we can 

show that that they can be offset by a suitable rotation of the axis of the 
symmetry I C C. Indeed, rotation of I by an angle a results in the rotation 
of S\ by the angle 2a. Yet for any two m-point sets A, B C {t = 1} on 
the unit circle, one can always find a rotation p of the circle so that the 
arclength distance between the points of A and p(B) will not exceed 2ir/m 2 
(a simple consequence of the pigeonhole principle). This shows that even 
in the worst case when the singular locus belongs to a single circle, the 
distance between S and can always be assumed bounded from below. 
The general case of arbitrary singular locus can be immediately reduced 
to this case by controlling the angular distance between the loci and their 
distance from the center. 

4.5. Demonstration of Theorem 2 

Consider an arbitrary Fuchsian linear system (1.1) from the Fuchsian 
class J^ n ,m with the matrix form fi the corresponding system of linear 
equations in the affine chart t as in (3.5). 
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Figure 4.1. "Thick slits" and genuine slits. 

This system after reduction to a linear differential equation of order n 
takes the form (3.6) with the coefficients a ,...,a„ G C[t] of the corre- 
sponding differential operator D of degree d < n 2 m and the slope not 
exceeding i? 2 ° ( " m) , R — R\,, by Lemma 4. Solutions of this equation are 
analytic outside the singular locus E = {t , . . . , r m }, though the leading 
coefficient a vanishes at some other points as well (apparent singularities 
of the equation) . 

Moreover, for any singular point r G S, one can find a real line r G 
t C C ("axis of symmetry") passing through this point, such that the 
linear operator D 6 obtained by symmetrization of D in this axis has the 
order 2n and the slope bounded by the same double exponential expression 
R 2 ' ' . Indeed, one has to apply Lemma 4 to the symmetrized system 
(4.15), in which fit is the result of reflection of fi in I and note that the 
carpeting function will be replaced by a polynomially equivalent one (4.16). 
Abusing the language, we will say that the operator D admits an axis of 
symmetry. 

Consider the system of "thick slits" of the disk Kr = {\t\ < R} as follows 
(see Fig. 4.5): 

(1) annuli {j ^ \t — Tj\ ^ \} around each singular point Tj G S; 

(2) rectangular strips of width | along all segments [t^Tj] (these seg- 
ments should be axes of symmetry for these rectangles); 

(3) rectangular strips of with | along the segment connecting tj with 
the exterior circle 8Kr = {\t\ = R}. 
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By elementary geometric! considerations it is clear that one can construct 
a system of genuine slits of K R \ S as follows: 

(1) these slits are only circular arcs and line segments; 

(2) they subdivide Kr \ S into simply connected polygonal domains 
and punctured disks; 

(3) each slit (circular or straight) is at least r-distant from the null locus 
of the leading coefficient V = {a = 0}, where 1/r = 0(n 2 m); 

(4) the total length of the slits is bounded by 0(m 2 R). 

Indeed, even if all d — m 2 n points from V belong to a given annulus or 
a rectangular strip of width O(l), one can always find a circle (resp., line 
segment) which is at least 0(l/d)-distant from V. 

We can now assemble all bounds given by Lemma 4, Lemmas 8 in com- 
mon terms relevant to the Fuchsian class ^ n ^ m and the carpeting func- 
tion R = R\, on it. The variation of arguments of any linear combination 
of solutions of the system is explicitly bounded along all slits constructed 
above, and using the argument principle for simply connected pieces (resp., 
Lemma 10 for annuli), we arrive at the explicit bounds for the counting 
function. 

Obviously, in all products (4.3), (4.6) and hence in (4.10), the double 
exponential term S = R 2 ° {d> absorbs all other terms without changing the 
asymptotic behavior of the entire products, which ultimately proves the 
bound asserted in Theorem 2. □ 
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